Theory of NMR as a local probe for the electronic structure in the mixed state of the 

high-T c cuprates 
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We argue that nuclear magnetic resonance experiments are a site-sensitive probe for the electronic 
spectrum in the mixed state of the high-T c cuprates. Within a spin-fermion model, we show that 
the Doppler-shifted electronic spectrum arising from the circulating supercurrent changes the low- 
frequency behavior of the imaginary part of the spin-susceptibility. For a hexagonal vortex lattice, 
we predict that these changes lead to (a) a unique dependence of the 63 Cu spin lattice relaxation rate, 
l/Ti, on resonance frequency, and (b) a temperature dependence of Ti which varies with frequency. 
We propose a nuclear quadrupole experiment to study the effects of a uniform supercurrent on the 
electronic structure and predict that Ti varies with the direction of the supercurrent. 



PACS numbers: 74.25.Nf, 74.60.Ec, 74.25. Ha 

The form of the fcrmionic excitation spectrum in and 
around a vortex in the mixed state of high-T c cuprates 
has been the subject of intense research in the last few 
years |Q0|- It has been known for a long time that the 
supercurrent circulating a vortex gives rise to a Doppler- 
shift in the fermionic excitation spectrum [pj. For a d- 
wave order parameter, Volovik pointed out [HJ, this shift 
leads to the scaling relation D{ep) ~ vH, where D(cf) 
is the density of states (DOS) at the Fermi energy, and H 
is the applied field. The main contribution to the DOS 
comes from regions close to the nodes, but in real space 
far from the vortex cores. Specific heat (SH) measure- 
ments by Moler et al. Q on YBa2Cu3 06.95 confirmed 
this scaling behavior, however, these experiments only 
provide information on the sample averaged DOS. 

In this communication we argue that the spin-lattice 
relaxation time, Ti, measured in nuclear magnetic reso- 
nance (NMR) experiments, is a more direct probe since 
it can provide site-specific information on the electronic 
structure in the vicinity of a vortex. Our conclusions are 
based on two important results. First, it was recently 
pointed out by Curro and Slichter (CS) that the depen- 
dence of resonance frequency on the local field allows one 
to measure T\ as a function of the nucleus position in a 
vortex lattice In addition, T\ -measurements are a 
purely electronic probe, in contrast to SH experiments, 
which arc dominated by phononic contributions. Second, 
within a spin-fermion model Q , we find that the Doppler- 
shifted electronic spectrum significantly changes the low- 
frequency behavior of the imaginary part of the spin sus- 
ceptibility, x"(Qi w )- Since x"(Q: w ) changes with the 
spatially varying supercurrent, l/Ti which is a weighted 
momentum average of x" m the zero frequency limit, ac- 
quires a frequency, i.e., spatial dependence. 

Experimentally, CS found for the 17 resonance in 
YBa 2 Cu307 that l/Ti increases with decreasing distance 
from the vortex core ||], i.e., with increasing frequency. 
Preliminary results of Milling et al. |?J show a qualita- 
tively similar behavior for the 63 Cu relaxation rate. We 
demonstrate that the experimentally observed frequency 



dependence of T\ can be qualitatively explained by the 
spatial variation of the supercurrent and the resulting 
Doppler-shift in the electronic spectrum and predict a 
temperature dependence of the relaxation rate which 
changes uniquely with frequency. Both results provide 
insight into the local structure of electronic and magnetic 
excitations in the mixed state, in contrast to the sample 
averaged information provided by SH experiments. 

In what follows we consider for simplicity a single-layer 
system; an extension of our results to the case of a double- 
layer system, i.e., YBa2Cu30g+a:, is straightforward. The 
63 Cu spin-lattice relaxation rate, 1/T±, for an applied 
field parallel to the c— axis is given by 
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where 



F c (q) = A ab + 2B(cos(q x ) + cos(q y ] 



(1) 



(2) 



and A a b and B are the on-site and transferred hyperfine 
coupling constants, respectively ||||. To describe the 
effects of a supercurrent on x" we use the spin-fermion 
model H in which the damping of the spin mode is de- 
termined by its coupling to the planar quasi-particlcs. 
Within this model, the spin propagator, x> is given by 
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where xo is the bare propagator, and n is the bosonic 
self-energy given by the irreducible particle-hole bubble. 
It was earlier suggested [0 that Re% -1 should possess 
a hydrodynamic form such that for lu — ► 
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where £af is the magnetic correlation length, a is a tem- 
perature independent overall constant, and Q = (tv, tt) 
is the position of the magnetic peak in momentum space 



1 



which we assume to be commensurate pj . The descrip- 
tion of Re hi Eq.(||) as a relaxational, non-dispersing 
spin mode at low frequencies is in agreement with the 
analysis of NMR and inelastic neutron scattering ex- 
periments on several I^-xSr^CuC^ and YBa2Cu306+ x 
compounds ||^ Jl0l , [l2"f . 

We are thus left with the calculation of the imagi- 
nary part of II which describes the spin damping brought 
about by the decay of a spin excitation into a particle- 
hole pair. In the superconducting state we find to lowest 
order in the spin-fcrmion coupling g 
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where G and F are the normal and anomalous Green's 
functions. Assuming that the supercurrent momentum, 
p s , possesses only a weak spatial dependence, we find in 
semiclassical approximation B 



G(k, in r , 



F(k, i£l m ) = -u k v k < - 



1 



1 



iCl m — Ek ifl m + i?_k 



(0) 



Here, E k is the dispersion of the fermionic quasi-particles 
Bogoliubons), which up to linear order in p s is given by 



E k = 



|A k | 2 +v F (k)-p a 



(7) 



where ek is the electronic normal state dispersion and 
v^(k) = de^/dk. Any changes in the <i-wave gap Ak = 
Ao (cos(k x ) — cos(fcj,))/2 and the BCS coherence factors 
Wk, Vk due to a supercurrent appear only to order p 2 and 
can thus be neglected. 

Since the dominant contribution to Imll for oj — > 
comes from quasi-particle excitations in the vicinity of 
the nodes, we expand E k around the nodes to linear order 
in momentum and perform the momentum and frequency 
integrations in Eq. (|5|) analytically. Combining the results 
with Eqs. ([!]), (|J) and (^) we find that the dependence of 
1/Ti on temperature T and on p s is determined by the 
set {D m /T}, where D m = Vp p s and Vp is the Fermi 
velocity at the node in the m'th quadrant of the Brillouin 
zone. The full expression for 1/Ti for arbitrary {D m /T} 
is too cumbersome to be presented here, however, in the 
low temperature limit where \D m /T\ 3> 1 for m = 1..4, 
it simplifies and we obtain up to order (T/D m ) 
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where the sum runs only over those nodes with D m , D n < 
0, C = {ag) 2 k B {h 2 lnle ) 2 /(2K) 1 v A = \8A k /dk\ at the 
nodes, and 
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Here, q n . m is the momentum connecting nodes n and m. 
The opposite signs of the T 2 -terms on the r.h.s. of Eq.(||) 
arises from the combination of Fermi functions involved 
in different relaxation processes. While the first term 
describes a process in which a Bogoliubon is simply scat- 
tered at the spin-fermion vertex, the second one involves 
processes in which two Bogoliubons are simultaneously 
created or destroyed at the vertex. Since the dominant 
contribution to the relaxation rate in the low tempera- 
ture regime comes from the second term in Eq.(^|), we 
find the unexpected result that in this limit \jT\T ac- 
tually decreases with increasing temperature. Note, that 
for T -> it follows from Eq.(|) that 1/TiT - p 2 s . 

In the high-temperature limit (\D m /T\ <C 1 for all m), 
we obtain to leading order in (D m /T) 2 
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where the sum runs over all nodes. For T 3> |T> m | we 
thus recover as expected the temperature dependence of 
the relaxation rate in the absence of a supercurrent. 

We now turn to NMR experiments in the mixed state 
of the high-T c cuprates which possess a hexagonal vortex 
lattice. We consider for definiteness YBa2Cu3C>7, where 
v F « 0.4 eV, v A ~ 20 meV, A ab /B = 0.7, and Uf ~ 2 
is temperature independent in the superconducting state 
Q . Each nucleus in the sample is characterized by a res- 
onance frequency Az/(r) = 7„?iiJ z (r) and a supercurrent 
momentum p s (r) . Here H z (r) is the local magnetic field 
and 7 n is the nuclear gyromagnetic ratio. Since the elec- 
tronic Zeeman-splitting for typical applied fields is of the 
order 10 _1 meV, and thus smaller than the Doppler-shift 
for most of the nuclei, it can be neglected. 

Our scenario for the calculation of T\ is only applicable 
to a given nucleus, if in its vicinity A(r) is uniform and 
p s (r) varies sufficiently slowly. Since both assumptions 
are violated in the vicinity of a vortex core we exclude in 
what follows all nuclei within a radius R — 2£ £, « 6ao 
of the vortex core. Here, £ a b and ao are the supercon- 
ducting in-plane coherence length and the crystal lattice 
constant, respectively, and H z (r) is given by 
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where A is the magnetic penetration depth, <E>o is the flux 
quantum, and q runs over the reciprocal lattice of the 2D 
vortex array. Since non-local as well as non-linear effects 
can effectively be accounted for by a renormalization of A 



2 



63 Cu $ 



T=1 .OK § 




48.95 48.975 49.00 49.025 49.05 49.075 49.10 49.125 

H (in units of tyj'fi) 

FIG. 1. The 63 Cu spectrum (solid line) and 1/TiT (open 
squares) as a function of the local field for T=l K. Both curves 
are truncated at a distance 2£ a b from the vortex core. 

jl3| we neglect them in the following. We can then calcu- 
late the supercurrent momentum via V x H(r) ~ p s ( r ) 5 
and thus obtain a spatial relation between H z (r) and 
p s (r). In Fig. [l] we plot the resulting 63 Cu spectrum 
(solid line) and 1/TiT (open squares) as a function of the 
local field H z at T — 1 K. Nuclei at the highest frequen- 
cies are located at a distance 2£ a {, from the center of the 
vortex, nuclei at the lowest frequencies are in the center of 
a vortex triangle, those at the maximum of the spectrum 
are at the midpoint between two vortices. For T = 1 
K, the relaxation rate for practically all nuclei is deter- 
mined by the low-temperature expression, Eq.(^j). In this 
case, 1/TiT ~ p 2 and hence reflects the frequency depen- 
dence of p s which increases with decreasing distance from 
the vortex core. The Doppler-shifted fcrmionic spectrum 
thus gives rise to a unique dependence of 1/TiT on fre- 
quency within the vortex-lattice spectrum of 63 Cu. This 
result is in qualitative agreement with preliminary mea- 
surements by Milling et al. jjj who indeed find that l/TjT 
increases with increasing local field. 

Another unique signature of Doppler-shifted excita- 
tions appears in the temperature dependence of the re- 
laxation rate. In Fig. || we present 1/TiT as a function 
of field for three different temperatures (the 30 K and 60 
K curves are horizontally offset for clarity) . For nuclei at 
the low-field end of the curve, p s is small, \D m /T\ <C 1 
already at low T, and hence (TiT) -1 ~ T 2 which in- 
creases with increasing temperature. On the other hand, 
for nuclei closer to the vortex core, and particularly for 
those at the high-field end of the curve, p s is large and 
\D rn /T\ > 1 even for T = 60 K (which is still below 
the lattice melting temperature at T m = 70 K). The 
temperature dependence of 1/TiT is therefore described 
by Eq.(||), and thus decreases with increasing tempera- 
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FIG. 2. (1/TiT) as a function of the local field for three 
different temperatures. The 30 K and 60 K curves are hori- 
zontally offset for clarity. 

ture, as discussed earlier. Between these two extrema 
there exists a cross-over region, characterized by a min- 
imum in the relaxation rate. This non-monotonic be- 
havior of l/TjT arises from two competing contribu- 
tions: while quasi-particle excitations involving nodes 
with D m > become exponentially suppressed with in- 
creasing H z , the contributions from the remaining (un- 
suppressed) excitations scale as ~ pi and thus increase. 
The cross-over occurs approximately at those nuclei for 
which max\D m \ ~ 0(T)\ the minimum thus shifts to- 
wards higher fields, i.e., larger values of |TJ TO |, with in- 
creasing temperature, as shown in Fig. 0. 

While the absolute scale of our results in Figs, [l] and 
H depends on a number of parameters, it is practically 
independent of ^af for £af > 3. In contrast, the pre- 
dicted temperature and frequency dependence of Ti de- 
pends only on {D m /T} and thus presents a generic fea- 
ture of a Doppler-shifted electronic spectrum due to the 
distribution of supercurrents in a hexagonal vortex lat- 
tice. 

An alternative experiment which could observe the ef- 
fects of a Doppler-shifted fcrmionic dispersion is a nuclear 
quadrupole resonance (NQR) measurement in which a 
uniform supercurrent is applied to the sample. Since 
in the low-temperature limit, Eq.(||), \jT{T depends on 
D m = Vp ■ p Sl with Vp being fixed by the under- 
lying lattice, one would expect that the relaxation rate 
varies with the angle between the supercurrent momen- 
tum, p s , and the crystal axes. In Fig. || we plot 1/T]T, 
for |p s | = 0.02 as a function of the angle <j> between p s 
and the crystal x— axis (see inset). For <f> = 0, i.e., for 
p s parallel to the a;— axis, we find T>2,3 = —vp q s /^/2 
while D 1)4 = +v F q s j\[2. 1/TiT thus decreases with T 
for \D rn /T\ > 1, and increases as ~ T 2 for \D m /T\ < 1, 
in agreement with the results in Fig. |j| On the other 
hand, <p = 7r/4, i.e., p s along the lattice diagonal, is a 
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FIG. 3. \jT\T for p s = 0.02 and three different tempera- 
tures as a function of the angle cf> between p s and the crystal 
a;— axis (see inset). 



special case since D1.3 = ±vf<Is while Z?2,4 = 0. While 
the contribution to 1/T\T from quasi-particle excitations 
involving only nodes 2 and 4 scales as T 2 for all temper- 
atures, the contributions involving only nodes 1 and 3 
decreases with increasing T. As a result, \jT{T increases 
for all temperatures, though only weakly at low T. For 
T 3> vp q s , the relaxation rate for all (f> scales as T 2 and 
thus becomes angle- independent. Note, that the strong 
angular dependence of the relaxation rate for intermedi- 
ate temperatures explains the spread in values of 1/T\T 
at a given local field which we found in Fig. |^, particu- 
larly in the T=60K curve. 

Spin-diffusion and vortex vibrations can poten- 
tially contribute to the spin-lattice relaxation and thus 
smear out the effects described above. However, spin- 
diffusion is strongly suppressed by the inhomogeneity 
of the magnetic field in a vortex lattice ]15| ], More- 
over, calculations on the effect of vortex vibrations in 
YBa2Cu306+:c p5[ as well as an experimental compar- 
ison of 17 and "Cu relaxation rates || conclude that 
vortex vibrations are irrelevant in the relaxation of 63 Cu 
but may play an important role in that of 17 0. 

Recently, several groups have calculated T± in the su- 
perconducting state starting from a Fermi gas (FG) de- 
scription of the cuprates, i.e., neglecting antiferromag- 
netic spin fluctuations |l5|-|l7|]. We extended our calcu- 
lations to the FG limit by using Xfuii — n/g 2 and found 
two distinct differences between the relaxation rate for a 
FG and that in the presence of strong spin fluctuations. 
First, 1/T\T for a FG always increases with temperature 
as long as A a b and -Bpossess the same sign, in contrast 
to our results in Eq.(g). Second, in the FG limit \jT{F 
increases monotonically with increasing local field in a 
vortex lattice, i.e., it does not exhibit the local minimum 
we found in Fig. |L We propose that these differences in 
the predicted relaxation rates enable NMR experiments 



to determine which of the two limits applies to the su- 
perconducting state of the high-T c cuprates. 

Finally, we found above that the behavior of the relax- 
ation rate in a vortex lattice reflects the presence of a su- 
percurrent, magnetic fluctuations and nodes in the super- 
conducting gap. We are currently studying |l8| whether 
similar effects also occur in Sr2Ru04 for which strong 
indications of ferromagnetic fluctuations and a p-wave 
order parameter exist p9| [. 

In summary, we have demonstrated that NMR is a site- 
specific probe for the electronic structure in the mixed 
state of the high-T c cuprates. We have shown that in a 
hexagonal vortex lattice, the Doppler-shifted electronic 
spectrum gives rise to a characteristic temperature and 
field dependence of the 63 Cu relaxation rate. We propose 
an NQR experiment in which the direction of a uniform 
supercurrent with respect to the crystal lattice is varied, 
and predict a unique angular dependence of T\ . Finally, 
we argue that our strong coupling results are qualita- 
tively different from those predicted for a FG. 
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